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ABSTRACT 

Approximation formulas are found for x(0) and x(l 
x(t) satisfies x = f(x, t) , x(0) = x Q , x(l) = x 1# The 
are applied to an example of two-body motion. 
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SOLUTION OF LAMBERT'S PROBLEM 
FOR SHORT ARCS 


1. INTRODUCTION 

Consider the following boundary- value problem: 


X = f(x, t) , x(0) “ x 0 , x(l) - X t , (1) 


where t is a scalar, x and f(x, t) are column matrices, and the over-dots in- 
dicate differentiation with respect to t. 

Formulas will be found for x 0 and x l such that 

* 0 = <£(°) > = ^C 1 ) . 

where $(t) satisfies 


o 

II 

X 

o 

0(0) 

“ x(0) , 

n 

/-n 

O 

x(0) 

0 ( 1 ) = x(l) , 

&1) 

= *(1> > 

11 

x(l) 


It will not be necessary to find 0(t) but only to make certain assumptions as 
to its form. 

The problem defined by (1) has an interesting history in celestial mechanics, 
going back to Euler, Lambert, Lagrange, and Gauss. Recent books by Battin [1] 
and Escobal [2] describe a total of eight methods for its solution in the case of 
an inverse- square, central force field. All the methods, however, are iterative, 
requiring considerable computation to obtain a solution. The method developed 
in this paper, while limited to cases of moderate time-span, offers a concise, 
explicit formula with no need for iteration. Within this time- span limitation, 
the method has the further advantage of allowing for completely general velocity- 
independent force functions. 
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The author has found the method useful in the calculation of short arc in- 
tercept trajectories for space vehicles, in preliminary orbit determination from 
observations of the position of a spacecraft at two times, and for ephemeris 
interpolation. 


2. AN APPROXIMATE SOLUTION 


Let a be any component of x, a(0) 
mations of the form 


a o> a (l) " a i> and assume approxi- 


“o " a 0 a 0 + b 0 a i * a 2^0 b 2*i + a 3* (i ’o + b 3<\ - 


( 2 ) 


a x - c 0 a 0 + d^ + c 2 a 0 + d 2 a 1 + c 3 a 0 + d 3 'a' 1 . 


( 3 ) 


The scalars a 0 , b 0 , a 2 , b 2 , a 3 , and b 3 are determined by assuming (2) to be 
exact when a ~ (t),i = 0, . . . , 5, where the <jL's are linearly independent 

over the interval [0,1] with derivatives through the third order at t = 0 and 
t = 1. The coefficients in (3) are determined in a similar way. If we use the 
same set of s for each component of x, we can write 


x, 


a 0 X 0 + b 0 X l + a '2 X 0 + b 2 x i + a 3 " X 0 + b 3 x i > 


( 4 ) 


c o x o + d o x i + c 2 x 0 + d 2 x l + c 3 x 0 + d 3 x 1 . 


( 5 ) 


We eliminate the third derivatives in (4) and (5) by using 


• •• _ * 

x - Px r w , 


( 6 ) 


obtained from (1), where P is a matrix with element in the ith row and jth column 
equal to the value of <9 f L /dx 3 , and w is a column matrix with ith element equal to 
the value of df l /dt, f 1 and x J being respectively the ith component of f and the 
jth component of x. Substituting (6) into (4) and (5), we obtain 

a 3 ^o) X 0 b 3^1 X l ” ft > (7) 
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C 3^0*0 + ^3^l)*l 


” y , 


(8) 


where I is the unit matrix and 


(3 - a Q x 0 + b oXl + a 2 x 0 + b/x, + a 3 w Q + b 3 w t , 

y = c 0 X 0 + d 0 X l + c 2 *0 + d 2*l * C 3 W 0 * d 3 w ! ’ 

x 0 and x’j being computed from x = f(x, t). Solving (7) and (8) we obtain 


(B + bP,P 0 ) x„ 

= £ + P, (b 3 r-d 3 /3), 

(9) 

(B + b P„ P, ) i, 

= ’ ,tp o( c j( 3 " a ) 5')' 

(10) 


B = l - a 3 P 0 - d 3 P lt 

b = a 3 d 3 " b 3 C 3 * 


Note that B + b P 0 P x is the transpose of B + bPj P 0 if P is symmetric, as is the 
case for a central force field. 


3. A POLYNOMIAL APPROXIMATION 
Substituting successively 

a = <£. (t) = t 1 , i = 0, 1, 2, 3, 4, 5 

into Equation (2) , we obtain 

a 0 + b 0 = ° , b 0 = 1 , 
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2a 0 

+ 2b 2 E 

■ 0 

l 

6b 2 

+ 6a 3 + 

6b 3 

= 0 , 

12b 2 

+ 24b 3 

- 

0 , 

20b 2 

+ 60b 3 

= 

0 . 


The solution of this set of equations is 


1 i kg 1 i a 2 ~ 


1_ 

20 > b 2 = 


JL 

20 - a 


3 “ 20 ' b 3 ~ 30 * 


In a similar way we find 


- 1 . d = 1 f 


3 

20 


CL ~ nr\ 
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20 ' 


C 3 “ 30 ' d 3 


■20 -(12) 


4. NUMERICAL EXAMPLE 

For a particle moving in an inverse- square central force field, (1) becomes 


x - 


M(r,-r 0 ) 


x ( 0 ) = x Q , 


x(l) = Xj , (13) 


where x is the position vector of the particle, fi is a constant, the real time r and 
the pseudo-time t being related by =r-r 0 with x(r Q ) = x Q , x(Tj) = Xl . 

As a numerical example take two vectors from a circular orbit separated by 
an angle of 15°. We can let r Q = o, r x = 77/12 fja, and 



T 


cos 15 ° 

x o " 

_o_ 

- x i “ 

sin 15 ° 
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Substituting these values and (11), (12), (13) into (9) and (10) gives 


-.00000022 


-.06775845 

261799360^ 

, x, = 

_. 2528788 19_ 


The true values are 



0 


-.06775867 


^261799388_ 

- x, - 

. 252878790 


5, REMARKS 

The method developed above can be carried through when (1) has the more 
general form 


x - f ( x , t) + A(t) x , x(0) - x Q , x(l) = x t , 


where A( t ) is a matrix function of t . 

The time-span restriction may be less severe when x in (1) represents a 
small deviation from a reference trajectory, where the solution of the boundary 
problem for the reference trajectory has been found by other means. For such 
problems we have x Q = x l = 0. For example, we might obtain a reference 
trajectory by solving a two-body problem with the given boundary conditions by 
one of the standard methods [1,2]. Then we could apply the method developed 
above to Encke's formulation [1] of the differential equation for the deviation of 
the position vector from the reference trajectory. 

If we attempt to approximate x 0 and x l in such a way that x 0 = 0(0) and 
Xj = 0(1) , where 0(t) satisfies 0 (4) (0) = x (4) (0) and 0 (4) (1) = x< 4 ) (1) as 
well as the other relations specified in the introduction, we are forced to solve 
a set of zionlinear equations unless (1) is linear in x. 
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